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Abstract 

The  scattering  problem  for  an  inhomogeneous  two-dimensional  biisotropic  cylinder  is  solved 
in  the  frequency-domain  by  means  of  a hybrid  method,  in  which  finite  difference  equations 
in  the  interior  region  are  combined  with  a mesh  truncation  in  terms  of  a boundary  integral 
equation.  Numerical  results  for  the  bistatic  echo  widths  are  presented  and  compared  with  a 
reference  solutions  in  the  circular  cylinder  cases  and  it  is  found  that  the  method  yields  more 
accurate  results  than  what  can  be  achieved  with  a local  absorbing  boundary  condition. 

1.  Introduction 

A scattering  problem  involving  inhomogeneous  media  can  be  solved  using  a volume  integral 
equation  (VTE)  that  is  discretized  into  a matrix  equation  by  using  the  method  of  moments 
(MoM).  However,  a drawback  with  VIE  is  that  the  integral  operator  discretizes  into  a filled 
matrix,  that  becomes  very  large  and  exceedingly  time  consuming  to  invert  if  the  geometry  is 
several  wavelengths  in  size.  An  alternative  is  to  use  a partial  differential  equation  (PDE)  method, 
in  which  the  differential  operator  discretizes  into  a sparse  matrix  by  using  finite  differences  (FD) 
or  by  using  the  finite  element  method  (FEM).  However,  when  using  FD  or  FEM  in  an  open  region 
a radiation  condition  must  be  supplemented  by  truncating  the  computational  domain  with  an 
absorbing  boundary  condition  (ABC),  which  can  be  either  an  approximate  (local)  ABC  or  an 
exact  (global)  ABC.  One  example  of  a local  ABC  is  the  perfectly  matched  layer  (PML),  which  has 
been  extended  recently  to  include  bianisotropic  media  [1].  A global  ABC  is  usually  formulated 
in  terms  of  an  integral  operator  on  the  boundary  of  the  computational  domain,  which  yields  that 
the  FD  or  FEM  equations  then  must  be  solved  in  conjunction  with  a boundary  integral  equation 
(BIE).  Examples  of  such  combined  methods,  referred  to  as  hybrid  methods,  are  e.g.  FEM-BIE 
[2]  and  FDTD-BIE  [3].  In  this  paper,  we  consider  a hybrid  FD  FrequencyDomain-BIE  approach 
to  a two-dimensional  scattering  problem  for  an  inhomogeneous  cylinder  of  a biisotropic  material 
Interior  to  the  region  containing  the  biisotropic  cylinder,  we  use  finite  differences  and  on  the 
boundary  of  the  region  we  use  a BIE,  which  reduces  to  a contour  integral  equation. 

2.  The  Scattering  Problem 

Consider  an  in  free  space  located  biisotropic  scatterer  (see  figure  1)  described  by  the  following 
constitutive  relations  in  the  frequency-domain: 


D (r)  = e0e  (r)  E (r)  + ^/£^^  (r)  H (r) , B (r)  = ir)  H (r)  + y^oC  (r)  E (r) , (1) 
where  the  parameters  e,  fi,  £ and  ( are  arbitrary  functions  of  the  position  variable  r. 
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Figure  1:  The  scattering  configuration. 

The  scattering  problem  is  to  determine  the  scattered  fields  Es , Hs  when  the  scatterer  is  excited 
by  known  incident  fields  E',H\  created  by  an  extraneous  time-harmonic  (exp  (jwt))  current 
density  J;  see  figure  1.  Substituting  the  constitutive  relations  (1)  into  the  Maxwell’s  equations, 
we  eventually  arrive  at  the  following  equations  for  the  scattered  fields: 

V x Es  + jfco  {m oH*  + CES)  = -jko  ((»  - 1)  r]0H ! + C^)  , (2) 

V x r ftH*  - jfco  ( eEs  + H°)  = j k0  ((e  - 1)  E{  + imH[)  . (3) 

3.  The  FDFD  Equations  in  the  Two-dimensional  Case 

The  2D  FDFD  equations  are  derived  in  a similar  way  as  in  [4],  The  region  containing  the  scatterer 
is  divided  into  square  cells  (see  figure  2),  in  which  the  constitutive  parameters  are  approximated 
as  homogeneous  and  equal  to  their  true  values  at  the  mid  points.  In  the  outermost  layer  of  cells, 
the  parameters  are  e = fj,  = 1,£  = £ = 0.  The  FDFD  approximation  results  in  a linear  system 
of  equations  for  the  everywhere  continuous  ^-components  of  the  fields  [5]: 

inside  ~b  outside  "b  -MEH-bT^inside  = inside)  (4) 

inside  ~b  outside  "b  ■^HE-^'^jnside  “^inside!  (5) 

where  E®  jnside,  nside  are  evaluated  at  the  interior  FDFD  points,  indicated  with  filled  circles 
in  figure  2,  and  E®  outside,  JT|  0Utside  are  evaluated  at  the  extruding  FDFD  points,  indicated 

with  blank  circles  in  figure  2.  Meh5  Mhe  and  MgH  are  sparsely  filled  square  matrices, 

describing  the  FDFD  operators  at  the  interior  points  and  the  matrix  B describes  the  FDFD 
operators  on  the  extruding  points.  Jj^side  and  J?side  are  excitation  vectors,  obtained  as  linear 
combinations  of  the  components  of  the  incident  fields  at  the  interior  points. 
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Figure  2:  One  corner  region  of  the  computational  domain.  T is  the  contour  at  which  the  global 
absorbing  boundary  condition  is  calculated. 
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4.  Absorbing  Boundary  Condition  in  Terms  of  a BIE 


With  all  their  sources  inside  the  contour  T,  depicted  with  a thick  solid  line  in  figure  2,  the 
^-component  of  the  scattered  fields  satisfy  the  integral  equation 


IfjM  + £ [a  (r y)  §5  M - k M §§  (r.O 


dT'  = 0,  r on  T, 


(6) 


where  Fz  denotes  or  iff  and  G(r,r')  — — |Hq^  (£o  |r  — r'|)  is  the  2D  free  space  Green’s 
function.  After  a simple  MoM  approximation  using  pulse  functions  and  point-matching  at  the 
points  indicated  with  crosses  in  figure  2,  we  obtain  the  following  relation  (the  numerical  ABC) 
between  the  scattered  field  and  its  normal  derivative: 


P FI  = Q dj*. 


(7) 


where  Ff  and  <9nFf  are  vectors  containing  the  field  and  its  normal  derivative,  respectively,  and 
where  P and  Q are  matrices.  To  connect  the  ABC  to  the  FDFD  scheme,  the  elements  in  F® 
and  <9nFJ  are  approximated  from  the  two  nearest  FDFD  points  on  each  side  of  F by  linear 
interpolations  and  central  differences,  respectively,  i.e. 


Ft 


F*  . _ 4-  Fs  • . , 
z,  outside  1 a, inside 


F3  . , — P5. 

S,  rps  ^ outside  inside 
n-F z ~ 7 j 


(8) 


where  h (the  mesh  parameter)  is  the  side-length  of  a cell.  With  (8)  in  (7),  we  obtain  a matrix 
relation  between  the  field  values  at  the  different  side  of  F: 


F)FZt  outside  insider  f^7fZ)0utsi(je  C Hz  , insider  0) 

where  the  square  matrix  D appears  to  be  well-conditioned  and  hence  invertible.  Using  (9)  in 
(4)  and  (5),  it  thus  follows  that  the  fields  at  the  interior  points  are  obtained  from  the  following 
matrix  equation: 

Mee  Meh 
Mhe  Mhh 

where 


tps 

z,  inside 

" 7® 

"inside 
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2, inside 
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inside 

(10) 


MEe  = M/ee  + BD-1C,  MHH  = MiIH  + BD-1C.  (11) 

Note  that  since  neither  B nor  C nor  D depend  on  the  properties  of  the  scatterer,  the  operation 
BD_1C  in  (11)  has  to  be  carried  out  only  once.  Thus,  for  a given  mesh  and  frequency  the  ABC 
can  be  precalculated  and  stored  for  further  usage  when  considering  different  properties  of  the 
scatterer,  as  long  as  the  scatterer  is  confined  within  the  boundary  T. 


5.  Numerical  Example 

In  this  one  example,  consider  a circular  chiral  cylinder  with  the  radius  = A,  where  A is  the 
wavelength  in  free  space.  The  constitutive  parameters  are  chosen  to  be  e = 2.6  — 0.2j,  //  = 
1.5  — 0.1j,£  = — C = — j«  = — j (0.6  — O.lj).  For  the  computational  region  we  consider  two 
different  mesh  sizes,  41  x 41  FDFD  points  with  spacing  h = 0.05A  and  81  x 81  FDFD  points 
with  spacing  h = 0.025A.  The  numerical  results  are  compared  with  a reference  solution  obtained 
by  means  of  an  expansion  in  cylindrical  eigenwaves.  The  cross-polarized  bistatic  echo  width 
Across  = ^ctoss  for  waves  impinging  in  the  direction  <j>  = 0 is  depicted  in  figure  3,  with  the 
solid  lines  for  the  reference  solution,  the  dotted  lines  for  h = 0.05A  and  the  dashed  lines  for 
h — 0.025A.  The  numerical  results  indicate  convergence  and  are  in  good  agreements  with  the 
reference  solution. 
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Figure  3:  The  cross-polarized  bistatic  echo  width  for  a chiral  circular  cylinder  with  the  radius 
= A and  the  parameters  e = 2.6  — 0.2j,  fj,  — 1.5  — O.lj  and  k = 0.6  — O.lj.  The  cylinder  is 
illuminated  in  the  direction  <f>  = 0°. 

6.  Conclusion 

The  scattering  problem  for  inhomogeneous  two-dimensional  biisotropic  cylinders  has  been  con- 
sidered. Finite  difference  equations  for  the  interior  region  as  well  as  a contour  integral  equation 
realizing  a global  absorbing  boundary  condition  have  been  derived  and  implemented  numerically. 

Computationally,  the  present  FDFD-BIE  method  is  much  faster  than  a moment  method 
using  a filled  matrix  of  the  same  size.  Hence,  the  present  method  can  be  used  for  solving  2D 
scattering  problems  involving  scatterers  of  intermediate  sizes  (i.e.  a couple  of  wavelengths)  in 
both  directions  whereas  MoM  under,  the  same  computation  time,  only  can  be  used  for  consid- 
erably smaller  or  thinner  structures. 

Since  our  global  ABC  is  unlikely  to  generate  artifacts  from  the  boundary  in  the  solution, 
the  usefulness  of  the  hybrid  method  as  a fast  solver  in  an  optimization  approach  for  solving  the 
corresponding  inverse  problem  is  of  interest  for  forthcoming  research. 
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